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We consider electron transport through a quantum dot described by the Kondo model in the 
regime of large transport voltage V in the presence of a magnetic field B with max(V,B) 3> Tk- 
The electric current I and the local magnetization M are found to be universal functions of V/Tk 
and B/Tk, where Tk is the equilibrium Kondo temperature. We present a generalization of the 
perturbative renormalization group to frequency dependent coupling functions, as necessitated by 
the structure of bare perturbation theory. We calculate I and M within a poor man's scaling 
approach and find excellent agreement with experiment. 
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The transport of electrons through quantum dots in 
the Coulomb blockade regime is strongly suppressed. In 
the case of an odd number of electrons on the dot, how- 
ever, when the dot carries a spin S, exchange coupling of 
S to the spin of the conduction electrons in the two leads 
(L, R) gives rise to a Kondo resonance at low tempera- 
ture T <C Tjc (Tk is the Kondo temperature). Electrons 
may then cross the dot via resonance tunneling, thus cir- 
cumventing the Coulomb blockade, as seen in a number 
of recent experiments [jj, |2| ■ While there is a wealth of 
methods like Bethe ansatz, conformal field theory, nu- 
merical renormalization group (RG) and self-consistent 
perturbation theory available to treat the Kondo model 
in equilibrium, most of these methods fail in the presence 
of a finite bias voltage V, for V > T K . 

The first theoretical treatment of this problem dates 
back to the 1960s 0] when the current through a Kondo 
impurity, in the presence of finite V and magnetic 
field B, was calculated in perturbation theory (PT) in- 
cluding leading logarithmic corrections. These early 
works also attempted a resummation of the logarithms 
within Nagaoka's or Abrikosov's approximations. All of 
these works, however, neglected salient non-equilibrium 
physics of this problem (see below). 

A crucial difference between a Kondo dot in the regime 
of large voltage V 3> Tk and the equilibrium Kondo 
problem is the presence of inelastic processes, associated 
with the finite current through the dot, down to the low- 
est temperatures 0, 0] . These processes destroy coher- 
ence on an energy scale T and prevent the full formation 
of the Kondo singlet resonance state. As discussed in 
for a conventional Kondo dot described by an An- 
derson model, T - V/\n{V/T K ) > T K for V > T K . 
In this regime, which we consider here, a perturbative 
treatment is sufficient to capture the developing Kondo 
correlations. As is well known from the usual Kondo 
effect, even in the perturbative regime it is necessary 
to resum the leading logarithmic terms in order to re- 
cover the scaling behavior. For the equilibrium case it 
is known how to achieve this in a systematic and con- 



trolled way by employing the perturbative renormaliza- 
tion group method. For the noncquilibrium Kondo prob- 
lem a comparable treatment has not been developed so 
far. A certain resummation of perturbation theory (PT) 
has been considered 0] within an Anderson model for the 
differential conductance G(V,B). However, the summa- 
tion of leading logarithms remained incomplete. 

In this Letter we present the first systematic study of 
renormalized PT for the nonequilibrium Kondo problem, 
considering in particular the influence of a magnetic field 
B. We focus on the regime V ^S> Tk, for any B. Our 
results are also applicable for large B 3> Tk, and arbi- 
trary V. In both cases the full formation of the Kondo 
resonance is inhibited, and a perturbative treatment in 
the exchange coupling J is possible. At first sight one 
might expect the nonequilibrium Kondo effect to be sim- 
ilar to the high temperature regime T Tk of the usual 
Kondo model. It differs, however, in three qualitative as- 
pects from the latter: (1) The occupation probabilities of 
the local spin are not thermal and have to be determined 
by solving a quantum Boltzmann equation. This leads to 
an unusual dependence e.g. of the spin susceptibility on 
V and to a novel structure of the logarithmic corrections 
in PT. (2) Electrons in a wide range of energies ~ V con- 
tribute to the low-T properties. In contrast to the case 
of large T, where all features are smeared on the scale 
T, we find that it is therefore essential to keep track of 
the frequency dependence of the renormalized coupling 
constants. (3) Life-time effects play a crucial role in non- 
equilibrium, as the corresponding rates are larger than T 
but smaller than V 

As a first step we present and discuss the results of 
bare PT for the current I(V, B) and the magnetization 
M(V,B) in leading logarithmic order. Next we turn to 
the derivation of renormalization group equations for the 
running coupling constants. The structure of bare PT 
forces one to keep the frequency dependence of the cou- 
pling functions. Employing simplifications in the spirit 
of Anderson's poor man's scaling, we derive a set of RG 
equations in one loop order. These have to be supple- 
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merited by a self-consistent determination of the relax- 
ation rate T, providing a cutoff of the RG flow. Finally, 
the physical quantities current I, magnetization M and 
differential conductance G are calculated to leading or- 
der in 1/ ln[max(V, B)/Tk] by inserting the renormalized 
frequency dependent coupling functions into the lowest 
order PT expressions. In this way, all physical quantities 
show scaling behavior in V/Tk and B/Tk, where Tk is 
the equilibrium Kondo temperature. 
We consider the Kondo Hamiltonian 

H = ^ ( e k - /ia)c^ kcr C Qkc r - BS Z 

+ 7} Ja'aS- (4, k , CT ,T CT vCak<r), (1) 

a,a ; — L, 7?, k,k' : a,cr' 

where ^l,r — ±V/2. S is the spin 1/2 on the dot 
and t are the Pauli matrices. We shall use the dimen- 
sionless coupling constants g d = NqJll = NqJrr and 
9lr = NqJlr — NqJrl (assuming a symmetric dot), 
where No is the local density of states. The local spin 
S = | Yl-yy f\ T i~i' f~t' i s represented by pseudo-fermions 
(PF) in the sector of the Hilbert space with ^ 7 /|/ 7 = 1. 

The electric current I = (e/h)JLRlva[D^ R (t,t)] 
through the quantum dot may be conveniently ex- 
pressed through the Keldysh component of the 
contour-ordered correlation function Dr,n(t,t') — 

(-i) 2 12w(Tc{c[«At)^cmAt)-m}) To ob- 
tain the leading logarithmic corrections, we evaluate 
the diagrams shown in Fig. which for T, V, B much 
smaller than the band cutoff D yield 

I = ^ E {^9lr) 2 [u 2 (V)+u 2 (V + B)+u 2 (V-B) 
-M(c(B)~c(-B))], (2) 

with u p (x) = x(l + 2pg d In j^J , (p = 1, 2), and c{B) = 

coth ^ [m {V + B) + 2Vg d In ^ + 2Bg d In ^ ] where 
all logarithms are implicitly cutoff by T . Appelbaum's 
expression for the current @ differs from © in that 
he assumed M to be given by its thermal equilibrium 
value. In the zero-field limit M vanishes and @ reduces 

t0 1 = ^k(i9lr) 2 SV (l + Ag d \n^j derived earlier in 
Refs. 0, @- To calculate the current for finite B, one 
needs to know the local magnetization M = nj — and 
at finite voltage this necessarily involves solving a quan- 
tum Boltzmann equation. In diagrammatic language 
this is nothing but the Dyson equation of the lesser PF 
Green function G<(uj) in steady state: G 7 (w)F 7 (lj) = 
Ti^{w)A 1 (u)) 1 where T 7 is the imaginary part of the self- 
energy and Aj the PF spectral function. The leading 
logarithmic corrections are obtained from the self-energy 
diagrams in Fig. yielding 




FIG. 1: Feynman diagrams for (a) Dlr, (b) PF self- 
energies and (c) vertices entering the 1-loop RG equation. PF 
(electron-) propagators are displayed as dashed (full) lines. 



coth^ [f + gl ui(S)] + g 2 LR (c(B) + c{-B)) 
T = 2g 2 + u 1 (B) + 2g 2 LR {u 1 (V + B)-u 1 (V-B)-2B} 

with g± = (g^±g 2 LR ). Note that the magnetization differs 
from the equilibrium result even to order g°, as has been 
derived independently in (see also Eq. JBJ) below) in 
disagreement with [UJ . Most notably, the magnetic sus- 
ceptibility changes from the usual l/T Curie law to a 1 /V 
behavior for V ^ T. Corrections of order g 3 InD to the 
collision integrals in the quantum Boltzmann equation 
result in corrections of order g InD to M. Thus, they are 
much larger than the usual g 2 In D corrections obtained 
in equilibrium. In the limit V — > 0, Eq. simplifies 
to the non- interacting result M — tanh(B/2T), as we 
have neglected the subleading g 2 In D-corrections which 
can be obtained by including the ReE shift of the PF en- 
ergy levels. The effects of ImE are more important and 
are discussed below. 

In the scaling regime, V,B <§; D, bare PT is not valid 
and has to be resummed even for small couplings g and 
in the weak-coupling regime V, B ^> Tk (see Fig. EJd) . 
The method of choice for such a resummation is the per- 
turbative RG, using the basic idea that a change of the 
cutoff D can be absorbed in a redefinition of the coupling 
constants g 0. 

A close inspection of Eqs. 12IMII indeed reveals that in 
the regime where, under renormalization, D gets smaller 
than V, it is not possible to absorb a change of D in a re- 
definition of the couplings. For example, logarithmic cor- 
rections to gLR in the denominator of Eq. J3J are propor- 
tional to 2 \n(D/V) for B — > while the analogous correc- 
tion in the numerator takes the form \n(D /V) + \n(D /T). 
This apparent breakdown of scaling is due to the fact that 
electrons in a finite energy window fi R < u> < /ii con- 
tribute to low-energy properties. How much their scatter- 
ing is renormalized will depend on their respective posi- 
tions within this window: upon renormalization the cou- 
pling "constants" acquire frequency dependences. Taking 
this properly into account, one recovers scaling. One way 
to derive perturbative RG equations in such a situation 
is to start from so-called "exact" RG equations for (in- 
dependent) one-particle irreducible Green functions by a 
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FIG. 2: Renormalized coupling constants g z , a (uj) (dashed 
lines) and g±(u) (solid line) for B = IOOTr- and V = 70T K . 
For these values Eqs. (00 yield T = 6A7T K . 

straightforward generalization of equilibrium RG meth- 
ods We will not follow this route here, but suggest 
a substantially simpler (but less systematic) approach. 
A considerably more involved real-time RG scheme has 
been proposed by Schoeller and Konig but has not 
yet been applied to the present problem. 

We start from the observation, that all logarithmic 
terms in next-to-leading order PT stem from the sim- 
ple vertex renormalizations shown in Fig. ^ when the 
real part « l/(ui ± B/2) of the PF Green function is 
convoluted with the Keldysh component of the electron 
line —2mN tanh[(oj — fj, a )/2T] . Using cutoffs symmetric 
with respect to (Al,r, respectively, one obtains at T = 



d 



dlnD 



sign u> 
jj — Au> 



2Q(D- |Aw|), 



(4) 



where Aw depends on fi a , B and the incoming and out- 
going frequencies. Now we include the effect of fre- 
quency and spin-dependent coupling constants denoted 
by g aCT " c ' Q J 7 " c for an incoming electron in lead a = L. R 

with energy uj c and spin a =1,1 interacting with a PF 
of spin 7 and frequency ojf describing the local spin. 
Primed quantities refer to outgoing particles. Generally, 
vertices are complex and have Keldysh indices, but we 
will only keep track of the real parts on one Keldysh 
contour, which we believe to be justified to leading order 
in 1/ \ii[V/Tk]- The pseudo fermion spectral functions 
are strongly peaked at u> = ±B/2, which allows setting 
ujf(id'f) to — 7-f ( — 7'"f )■ Furthermore, neglecting the re- 
maining frequency dependence of the coupling functions 
in the frequency window |cj| < D |l3|, one is led |b| to 
two RG equations for the coupling functions g za (uj) and 
g±(uj) (assuming J LL = J RR = J L r = Jrl = J) 



dlnD 

dg±(v) 
d\nD 



/3=-l,l 



7y ) V u+ r(B+0%) ( 5 ) 



- ( aB +f}V s~ (0V) 
9X\ 2 )9zcr{ 2 J q 



ui+- 



j8,or=— 1,1 



#7,-75/217,-7.8/2 



r9za(^) 



(6) 
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7,-7S/2;7,-7S/2 



(7^^ + 7^^)^(07-75/2), 



with the initial conditions g za (uj) — gj_(uj) — JNq at 
D = D , the bare cutoff, and 9 Aw = &{D - |Aw|). 
Despite the fact that we use the same functions for both 
diagonal (LL,RR) and off-diagonal (LR,RL) processes, 
their renormalization is drastically different as different 
frequency ranges are probed in the two cases (see below) . 
In various limits analytic solutions are possible but as 
the formulas turn out to be rather lengthy, we will here 
restrict ourselves mainly to numerical results. 

The RG equations (J5J are valid in the regime of cutoff 
values D where the effect of spin- flip processes destroying 
the quantum coherence of the Kondo bound-state is neg- 
ligible, or D > r. Due to the finite current, relaxation 
processes contributing to T are present even at T = 0, 
and lead to an imaginary part of the self-energy of the 
pseudofermions (e.g. in the PF propagator in the dia- 
grams of Fig. [IJ;) and to vertex corrections. Although 
slightly different relaxation rates may enter in the var- 
ious coupling constants, with logarithmic accuracy we 
may use a single rate T. As all relevant processes involve 
at least one spin-flip, we identify T with the transverse 
spin relaxation rate I/T2, which is given in terms of the 
renormalized couplings as 

r = h E WfciM 2 / u -„ji-/,-v) (7) 

„. _./_r o 



a, a — L,i? 
7=-l,l 



with g±(w) = g±(-u), g z \{u) = 9zl{-w) and 



+g ± (uj - 7 B/2) 2 f u ^ a (I - U-^,- 7 b) 



Note the missing factor 2 in front of the spin-flip term 
g±(u> — j B/2) 2 in T (compared to the corresponding term 
in jnj, below), which reflects the fact that I/T2 0] is not 
given by a pseudo-Fermion self-energy and arises from 
vertex corrections in a linear-response calculation of the 
transverse susceptibility. The effect of T2 is to provide an 
upper time cutoff for the coherent scattering processes 
entering Eq. JJJl. Hence, it is modeled phenomenologi- 
cally by replacing 6 W with Q(D — \J 1J 1 + T 2 ) in Eq. 

A self-consistent solution of Eqs. (0 01 can now be 
obtained with little numerical effort (formally, effects of 
self-consistency for T are subleading for V, B 3> Tk)- 
The resulting couplings g z n{u) are shown in Fig. [5] In 
the renormalization process the g z /±(w) develop peaks at 
frequencies ui — —o-B±^- for g ztJ and u> = ±i(l?±V) for 
g±, at which resonant scattering from one Fermi surface 
to another becomes possible. Only these resonant pro- 
cesses survive for small running cutoff D, as can be seen 
from the structure of the functions 8 W = Q(D — \u\) in 
ijSJl. This Kondo-type resonant enhancement is eventu- 
ally cutoff by r. Note that T remains larger than Tk for 
V,B> Tk prohibiting the flow towards strong coupling- 
contrary to the suggestion of ^(J (Ref- []| discusses condi- 
tions for strong-coupling physics in the regime V 3> Tk)- 
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FIG. 3: a) Local magnetization M(V/B, B/Tk) of a symmet- 
ric dot for fixed magnetic field B. b)Differential conductance 
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in units of Go = ^s. Comparison of perturbative RG and 
bare PT to 0(J 2 ) and 0(J 3 ). c) Conductance measurements 
of Ref. (symbols) on metallic point contacts in magnetic 
fields 0.85, 1.7, 2.55 T (B = 36, 72, 104 T K , with T K w 30mK 
0). Assuming that the corresponding point contact is de- 
scribed by a single-channel (Jlr = V JllJrr) Kondo model, 
^sa ~ 4.2 is determined from G(V = 0, 5 = 0, T = 50mK) = 
JJ^+jf^ G sy m(T /Tk), where G sy m is known exactly from 
NRG calculations [l6jl . This fixes a/Z parameters for our RG 
calculation (solid lines, T = 0) which uses a straightfor- 
ward generalization of for Jll Ji?.r. As the (in- 
dependent) background is not known experimentally, we sub- 
tract AG = G B - 5.2 • 10~ 5 Go^, where G B is fitted to our 
results at large V . The experimental temperature T = 50mK 
leads to an extra small broadening at V = B. 

Knowing the renormalized g z /±(uj) at D = 0, the mag- 
netization M can be found from solving the Boltzmann 
equation dtn a = or 



J—T. If* 



a,a'—L,R 



g±(w- B/2) 2 f^ a (l- U-^-b) (8) 



a,a'=L,R'' 



the solution of which in turn determines the current as 
I 



TTfi 

16ft 



/dw [feH'/^a-Z-pJ (9) 

« — i i 



<r=— 1,1 
\2 



+43 ± (w - aS/2) 2 /^- M i,(l - fu-^-aEiW 



In Eqs. 00 there are further lifetime effects, which be- 
come relevant for \V — B\, T < Y. This broadening we 
model phenomenologically by smearing the Fermi func- 
tions over a width V Note that expanding these 
results in bare PT, one recovers Eqs. (121: ill . 

In Fig.Elwe show M(V/B,B/T K ) and G(V/B, B/T K ). 
Within PT, G shows threshold behavior at \V\ = B, due 
to the opening of another transport channel involving 
spin flip, requiring the Zccman energy B. Higher order 
processes enhance the steps into peaks (Fig. |3Jd). The 
agreement with experiment is excellent (Fig. Eh), con- 
sidering that there is no free parameter and our result is 
valid up to terms of order 1/ ln[max(V, B)/Tk\ only. We 



emphasize that neither perturbation theory (Fig. ^p) nor 
Appelbaum's result (as shown in .2$), nor the result of [f| 
can describe the experiment. 

In conclusion, we propose a simple method to gen- 
eralize poor-man's scaling (one-loop order) to non- 
equilibrium. Compared to the equilibrium situation, it 
is necessary to include non-equilibrium distribution func- 
tions, decay-rates and the w-dependence of coupling con- 
stants. The different structures arising even in the per- 
turbative regime may be best exemplified by considering 
the local magnetic susceptibility for V ^> Tk, B — > 



X(V) = 



2a In y~ 



(1 



In 



V 



1 



^2a In ^)ln 



ln-v_(l +a l n .v_) 



(1+aln v ) ln v 



(10) 



with a = (g 2 d - g 2 LR )/{2g LR ) and T K = D e -V(»*+w«). 
To the same order of approximation, one obtains in equi- 
librium for T 3> Tk just x = l/T. 
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